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Abstract 

25th-order high-temperature series are computed for a general nearest-neigh- 
bor three-dimensional Ising model with arbitrary potential on the simple cu- 
bic lattice. In particular, we consider three improved potentials character- 
ized by suppressed leading scaling corrections. Critical exponents are ex- 
tracted from high-temperature series specialized to improved potentials, ob- 
taining 7 = 1.2373(2), V = 0.63012(16), a = 0.1096(5), r/ = 0.03639(15), 
(3 = 0.32653(10), 5 = 4.7893(8). Moreover, biased analyses of the 25th-order 
series of the standard Ising model provide the estimate A = 0.52(3) for the 
exponent associated with the leading scaling corrections. 

By the same technique, we study the small-magnetization expansion of the 
Helmholtz free energy. The results are then applied to the construction of 
parametric representations of the critical equation of state, using a systematic 
approach based on a global stationarity condition. Accurate estimates of several 
universal amplitude ratios are also presented. 
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1 Introduction 



The Ising model is one of the most studied systems in the theory of phase transitions, 
not only because it is the simplest nontrivial model that has a critical behavior with 
nonclassical exponents, but also because it describes the critical behavior of many 
physical systems. Indeed, many systems characterized by short-range interactions 
and a scalar order parameter undergo a continuous phase transition belonging to 
the Ising universality class. We mention the liquid-vapor transition in simple fluids 
and the critical transitions in multicomponent fluid mixtures, in uniaxial antiferro- 
magnetic materials, and in micellar systems. Continuous transitions belonging to the 
three-dimensional Ising universality class are also expected in high-energy physics, for 
instance in the electroweak theory at finite temperature and in the theory of strong 
interactions at finite temperature and finite baryon-number chemical potential. For 
a recent review, see, e.g., Ref. 

The high-temperature (HT) expansion is one of the most efficient approaches to 
the study of critical phenomena. Very precise results have been obtained by perform- 
ing careful extrapolations to the critical point, by using several different methods, 
see, e.g., Ref. 0. For moderately long series, such as those available for models in 
the three-dimensional Ising universality class, the nonanalytic confluent corrections 
are the main source of systematic errors. For instance, according to renormaliz at ion- 
group theory, the critical behavior of the magnetic susceptibility is given by the 
Wegner expansion 

^ = Cr^ (1 + a^t^ + a2 + ... + bt^' + ... + eit + ea + ...), (1.1) 

where t = {T — T^/Tc is the reduced temperature and A is a noninteger exponent, 
A 0.5 in the Ising case. In the analysis of HT expansions these nonanalytic terms 
introduce large and dangerously undetectable systematic deviations in the results. 

In order to obtain precise estimates of the critical parameters, the approximants 
of the HT series should properly allow for the confluent nonanalytic corrections . 
However, the extensive numerical work that has been done shows that in practice, 
with the series of moderate length that are available today, no unbiased analysis is 
able to take effectively into account nonanalytic correction-to-scaling terms. In order 
to treat them properly, one should use biased methods in which the presence of the 
leading nonanalytic term with exponent A is imposed (see, e.g., Refs. [p!0|-|r?|). An 
alternative approach to this problem consists in considering models — we call them 
improved — that do not couple the leading irrelevant operator that gives rise to the 
confluent correction of order t^. Therefore, such correction does not appear in the 
expansion of an?/ thermodynamic quantity near the critical point: for instance, a^^ = 
in Eq. (|1 . 1| ) . In this case, we expect standard analysis techniques to be much more 
effective, since the main source of systematic error should have been eliminated. 
There are no methods that allow to determine exactly improved models, and one must 
therefore use numerical techniques. One may use HT expansions, but in this case the 
improved model is determined with a relatively large error |jl|,p,||, |, |1^, so that 



the final results do not significantly improve the estimates obtained from standard 
analyses using biased approximants. Recently ||T9|-27|, it has been realized that Monte 
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Carlo (MC) simulations using finite-size scaling techniques are very effective for this 
purpose, obtaining accurate determinations of several improved models in the Ising, 
XY, and 0(3) universality classes. 

As shown in Refs. [ p!9| , [25| , p7| -p9[] , analyses of the HT series for the improved models 
lead to a significant improvement in the estimates of the critical exponents and of 
other infinite- volume HT quantities. Our working hypothesis is that, with the series 
of current length, the systematic errors, i.e., the systematic deviations that are not 
taken into account in the analysis, are largely due to the leading confluent correction, 
so that improved models give results with smaller and, more importantly, reliable 
error estimates. This hypothesis can be checked by comparing the results obtained 
using different improved models: if correct, they should agree within error bars. In 
the following we shall report results that confirm our hypothesis. Indeed, the esti- 
mates obtained from three different improved Hamiltonians are perfectly consistent. 
Moreover, they are very stable with respect to the order of the series considered in 
the analysis, without showing dangerous trends, but only an apparent reduction of 
the error. The results obtained in Ref. [^] using 20th-order series are fully consistent 
with the 25th-order analysis that we present. 

We consider a simple cubic lattice and scalar models with Hamiltonian 

n = -(3j2^^^^ + J2^(^')^ (1.2) 

<i,j> i 

where j3 = 1/T, <i,j> indicates nearest-neighbor sites, (pi are real variables, and 
V{(j)'^) is a generic potential satisfying appropriate stability constraints. These models 
are expected to have either a critical transition belonging to the Ising universality 
class or a first-order transition between a disordered and an ordered phase, apart 
from special cases that correspond to multicritical points. Using the linked- cluster 
expansion technique, we compute, for an arbitrary potential, the HT expansion of the 
two-point correlation function to 25th order on a simple cubic lattice. These results 



extend those of Ref. [^] that reported the two-point function to 20th order ||30[. In 
particular, we consider three classes of models depending on an irrelevant parameter, 
which is fixed by requiring the absence of the leading scaling correction. The first one 
is the 0^ lattice model with potential 

V{cf>') = <f + - If. (1.3) 

MC simulations using finite-size scaling techniques have shown that the model is 
improved for pi 

A4 = A4 = 1.10(2). (1.4) 



A consistent but less precise estimate can be obtained from the HT expansion p9| . 
The second class of models is the 0^ lattice model with potential 

V{<i?) = <i? + A4(0' - 1)' + A6(0' - If. (1.5) 

Fixing Ag = 1, the 0^ Hamiltonian is improved for [jT^ 

A4 = A4 = 1.90(4). (1.6) 
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Finally, we consider the spin-1 (or Blume-Capel) Hamiltonian 

n = -f3j2s.s, + Dj2sl (1.7) 

<i,j> i 

where the variables take the values 0, ±1. An improved spin-1 model is obtained 
for 

D = D* = 0M1{8). (1.8) 

The comparison of the results obtained using the above-mentioned improved Hamil- 
tonians represents a strong check of the expected reduction of systematic errors in 
the HT results, and provides an estimate of the residual errors due to the subleading 
confluent corrections to scaling. 

We also extend the HT expansion of the zero-momentum n-point correlation 
functions Xn- In particular, we compute X4; Xe, and xs to 21st, 19th, and 17th 
order respectively. The analysis of such series provides information on the small- 
magnetization expansion of the Helmholtz free energy in the HT phase. These results 
are used to determine approximate representations of the equation of state that are 
valid in the critical regime in the whole (t, H) plane. For this purpose, following 
Ref . |l^ , we use a systematic approximation scheme based on polynomial parametric 
representations and on a global stationarity condition. This approach allows us to 
obtain an accurate determination of the critical equation of state in the whole critical 
region up to the coexistence curve. 

In Table |I] we anticipate most of the results that we shall obtain in this paper. We 
report HT estimates of the critical exponents and of the coefficients parametrizing 
the small-magnetization expansion of the Helmholtz free energy: they are denoted by 
IHT, where the "I" stresses the fact that we are using improved models. Then, we 
report several amplitude ratios (definitions are given in Sec. Those appearing in 
the column IHT-PR are obtained from an approximate representation of the equa- 
tion of state that uses the HT results as inputs, those labelled by LT are obtained 
from the analysis of low-temperature expansions, while those reported under IHT- 
PR+LT are obtained combining the IHT-PR and LT results. The comparison with 
the corresponding Table XIII of Ref. [|T^] shows that the estimates obtained from the 
25th-order series are essentially identical to those obtained by using the shorter 20th- 
order series. However, the longer series allow us to give error bars that are smaller by a 
factor of 1.5-2, depending on the observable. The estimates reported in Table |l| are in 
substantial agreement with, and substantially more precise than, the best theoretical 
and experimental results that have been previously obtained p!0|,p0|-p^,[35|-|50|. For 
a comprehensive recent review of theoretical and experimental results, see Ref. [Q. 
On the experimental side, we mention the planned experiments in microgravity en- 



vironment described in Ref. which may substantially improve the experimental 
determinations of the critical quantities and make the comparison with the theoretical 
computations more stringent. 

After completion of this work, the study reported in Ref. appeared, where 
analyses of 25th-order series for spin-S* models are reported. Results for the critical 
exponents are obtained by means of biased analyses, essentially by fixing A. Com- 



paring Ref. |T3 with Refs. |jT0|,|13[, where 21st- and 23rd-order series are analyzed, a 
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trend appears towards better agreement with improved Hamiltonian results (Ref. ||T9[ 
and present paper). The latest results of the authors of Ref. are in full agreement 
with our estimates. 

The paper is organized as follows. 

In Sec. 1^ we report the HT expansions. Section |^ reports the results of our analysis 
of the HT series for the critical exponents. In Sec. | we determine approximate 
representations of the critical equation of state. In Sec. [4.1| we give the definitions, 
in Sec. |4.2| we give estimates of the zero-momentum four-point coupling and of the 
first few coefficients of the small-magnetization expansion of the equation of state, in 
Sec. [4.3| we explain the method, and in Sec. we give the final results. In Sec. |^ we 
present estimates of several universal amplitude ratios. In Sec. ^ we determine the 
low-momentum behavior of the two-point function in the HT phase. 



2 High-temperature expansion 

We considered a simple cubic lattice and computed the HT expansion of several quan- 
tities for a generic lattice model defined by the Hamiltonian ( |1.2| ), using the vertex- 



and edge-renormalized linked-cluster expansion technique, developed in Refs. [p!^,p2 



and described in detail in Ref. |5^. Some technical points that allowed us to ex- 
tend the computation of Ref. will be reported in a forthcoming publication. We 
computed the 25th-order HT expansion of the two-point function 

G{x) = (0(O)0(x)). (2.1) 

In the present context we consider its moments 

= ^|x|2^G'(x), (2.2) 

X 

and therefore, the magnetic susceptibility x = "^o and of the second-moment corre- 
lation length = 7712/ (6x)- 

We also calculated the HT expansion of the zero-momentum connected 2j-point 
correlation functions X2j 

X2j = (<^(O)<^(^2)--0(a;2i-i)0(a;2j))c (2.3) 

X2,...,X2j 

ix = X2)- More precisely, we computed X4 to 21st order, xe to 19th order, xs to 17th 



order. The correlation function Xio was computed to 15th order in Ref. |19 . 

It would be pointless to present here the full results for an arbitrary potential: 
the resulting expressions are only fit for further computer manipulation. They are 
available on request. In Table |^ we give the new coefficients only for the three 
improved models we have considered, i.e., for the 0^ model at A4 = 1.10, for the 



k6 



model at Ae = 1 and A4 = 1.90, and for the spin-1 model at D = 0.641. 
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Table 1: Summary of the results obtained in this paper (unless a reference is ex- 
plicitly cited) by our high-temperature calculations (IHT), by using the parametric 
representation of the equation of state (IHT-PR), by analyzing the low-temperature 
expansion (LT), and by combining the results of the two approaches (IHT-PR-I-LT). 
The estimates of critical exponents marked by an asterisk have been obtained using 



scaling and hyperscaling relations. 







IHT 


IHT-PR 


LT 


IHT-PR+LT 


critical 


7 


1.2373(2) 








exponents 


V 


0.63012(16) 










a 


0.110(2), *0. 1096(5) 










V 


*0.03639(15) 










P 


*0.32653(10) 










5 


*4.7893(8) 










A 


0.52(3) 










LO 


0.83(5) 












2.0208(12) 








small-magnetization 


at 


23.56(2) 








expansion of 


re 


2.056(5) 








the free-energy 


rs 


2.3(1) 








in the HT phase 


rw 


-13(4) |l|] 


-10.6(1.8) 






universal 


Uo 




0.532(3) 






amplitude 


U2 




4.76(2) 






ratios 






-9.0(2) 






see Sec. |5| 






0.0567(3) 






tor notations 






0.02242(12) 








D+ 

RJ 




7.81(2) 








V3 




6.050(13) 








i?4 




93.6(6) 








V4 




16.17(10) 








R 




1 660f4'l 








9 

W 






4.75(4) 














1.956(7) 








0.01880(8) 












0.2659(4) 








Q- 








0.00472(5) 




Qc 




0.3315(10) 








93 








13.19(6) 












76.8(8) 




Qt 


1.000200(3) 
















1.032(4) g 














1.896(10) 




Qt 








1.024(4) 




Q2 








1.195(10) 




p 




1.2498(6) 








Pc 




0.3933(7) 








Rp 




1.9665(10) 
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Table 2: Coefficients of the HT series for the improved models. Lower-order coeffi- 
cients appear in Ref. p9 |. 



A4 



1.10 



As = 1, A4 = 1.90 



spin-1: D = 0.641 



21 958465949.119795229380125 

22 2581828793.17418316658592 

23 6953921835.10625772660286 

24 18716342130.2600278822297 

25 50369768053.5367726030130 



X 

55356759.0258594943774739 
130996257.131383657648562 
309956395.981892002096689 
732873665.558914443007657 
1732674465.68758001711514 



521863527.549747127784405 
1367254366.70256684609648 
3581814299.63029965928082 
9376338630.49601545283933 
24543094928.9205155990856 



21 32990320251.5660972216018 

22 94071328367.8146359923071 

23 267461898855.689392585599 

24 758423675496.642760823002 

25 2145329356955.42924803892 



m2 

1900950559.23375555678011 
4762044317.91673448231502 
11894571003.1970044574018 
29631147101.2512233682029 
73634162230.2093808561076 



17908950773.4801706544197 
49684326561.5439542757331 
137433163639.457494472451 
379139772127.101469600055 
1043350926215.22611634874 



20 541141652908.631074719231 

21 1643345014677.80358819408 

22 4961021084766.33884428748 

23 14895796670810.3387628037 

24 44504475774409.2126174407 

25 132362288688779.709839376 



TO4 

35399348720.3598637148375 
94444621918.7858920241050 
250485298262.046958470064 
660748522303.208118944668 
1734347627024.93369651634 
4531641133142.45499870752 



299758906549.791610350073 
885976701269.736104292700 
2603026564263.78069815384 
7606210964865.32821158574 
22115153167519.1984380502 
64005596692608.8036008995 



19 -141558376231.985023846408 

20 -440895445559.088001425635 

21 -1363771989486.31756523825 



X4 

-9210343000.40488445467068 
-25206881115.0765162521666 
-68511054288.5805997438372 



-77210883309.3840433243811 
-234263398532.544236218037 
-705801443484.646787710146 



18 
19 


25922773662329.4681285982 
93214547843378. 1420243052 


X6 

1657400403425.39611029038 
5239283130720.37310719268 


13110582140461.8241625980 
46080008679021.7095625364 


17 


-3021378127745877.943411840 


-188904527250502.5683919596 


-1360671334948122.792253527 



For the standard Ising model, we give below the coefficients of the terms that 
extend the expansions presented in Refs. [0,^ for x, m2, and Xi- 

X = ■■■ + 18554916271112254i;2^ + 85923704942057238^2^ + 0(t;2^) 
m2 = ■■■ + 977496788431483776t;2^ + 4767378698515169334^25 + 0(^2^), 



= 6306916133817628^^^ - 34120335459595728t; 



19 

-183166058308506108^20 _ 976373577976196368i;2i + Oiv"^^), (2.4) 
where v = tanh (3. 
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3 The critical exponents 



In this section we shall report three different analyses for the critical exponents. In 



Sec. we shall use integral approximants and derive the estimates reported in Table 



|T|. In Sees. and ^]3| we shall use two other methods that have been recently used 



in the literature ||T3,|16|,|I^ to confirm the integral-approximant results. 



3.1 Analysis using integral approximants 

In order to estimate 7 and u, we analyze the HT series of the magnetic susceptibility 
and of the second-moment correlation length respectively. We follow closely App. B 



of Ref . |jT9[ , to which the reader is referred for more details. 

We use integral approximants (lA's) of first, second, and third order (see Ref. 
for a review). Given an nth-order series /(/?) = Y17=o'^iP^> A;th-order integral 
approximant [mk/rrik^i/ ■■■/mo/l] lAk is a solution of the inhomogeneous A;th-order 
linear differential equation 

Pk{f3)f^'\f3) + Pu^mf'-'\p) +■■■+ pmf'\P) + + R{(^) = 0, (3.1) 

where the functions Pi{(3) and R{(3) are polynomials of order rrii and / respectively, 
which are determined by the known nth-order small-/? expansion of f{f3)- Follow- 
ing Fisher and Chen we also consider integral approximants, FCIA/c's, in which 
Pk{(3) is a polynomial in FCIA/c's allow for the presence of the antiferromagnetic 
smg ularity at = -pc 0. In our analyses we consider diagonal or quasi- diagonal 
approximants, since they are expected to give the most accurate results. For each 
set of lAfc's we determine the average of the values corresponding to all nondefective 
lA/c's. The error bar from each class of lA's is essentially the spread of the results, 
and it is given by the standard deviation of the results obtained from all nondefective 
lA's. In most cases the nondefective lA's are more than 90%. 

All lA's considered give perfectly consistent results. Moreover, the results turn out 
to be very stable with respect to the number of terms of the series, so that there is no 
need to perform problematic extrapolations in the number of terms in order to obtain 
the final estimates. In Fig. |^ we show the estimates of 7 obtained by analyzing the 
series of x for the 0^ model at A4 = 1.10 by using lAl's, IA2's, IA3's, and FCIA2's, as 
a function of the order n of the series considered in the analysis. Perfect agreement 
is also found among the results for the three improved Hamiltonians. This is shown 
in Fig. ^ where the results of the IA2 analyses for the three improved Hamiltonians 
are reported versus n. In Fig. ^ we also show the results of the IA2 analysis applied 
to the series of x for the standard Ising spin-1/2 model. The corresponding results 
disagree with those obtained by using improved Hamiltonians: clearly, there is a large 
error that is not taken into account by the spread of the approximants. The results 
for the Ising model improve if one biases the analysis by using the very accurate MC 



estimate of /?c H: /?c = 0.22165459(10). Indeed, 7 drops from 1.245 to 7 = 1.2400(5). 
However, the error obtained from the spread of the approximants is still incorrect. 
Results that are closer to those obtained by using the improved Hamiltonians (and 
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1.248 



1.244 



1.240 



1.236 



1.232 



1 I I I I r 



I I r 

■ IA1's 

• IA2's 

^ FCIA2's 

▼ IA3's 




_L 



_L 



_L 



_L 



_L 



_L 



_L 



_L 



_L 



13 14 15 16 17 18 19 20 21 22 23 24 25 26 



Figure 1: Estimates of 7 as obtained by analyzing the HT series of x i^r the (p^ 
model at A4 = 1.10 versus the order n of the series considered in the analysis. Several 
approximants (defined in the text) are considered. 



substantially compatible with them) are only obtained by additionally biasing the 
series, allowing for 0{t^) confluent corrections, see, e.g., Ref. 
For the 0^ lattice model we obtained 

/?c(A4 = 1.10) = 0.3750975(5), (3.2) 
7e(A4) = 1.23732(10) + 0.006(A4 - 1.10), (3.3) 

where 7e(A4) is the effective critical exponent obtained in the lA analysis, which has 
a small but nonvanishing dependence on A4 around the favorite value A4 = 1.10. 
(Here and in the following, we report explicitly the dependance on A4 and equivalent 
couplings: should a better estimate of A4 become available, it can be immediately used 
to improve our results.) The number between parentheses is basically the spread of 
the approximants at A4 = 1.10. The A4-dependence is estimated by determining the 
variation of the results when changing A4 around A4 = 1.10. The best estimate of 7 
should be obtained at A4 = A4. Thus, using the MC estimate of A4, i.e., A4 = 1.10(2), 
and taking into account its uncertainty, we obtain the estimate 7 = 1.23732(10) [12] 
(which is also reported in Table ^), where the error in brackets is related to the 
uncertainty on A4. As final error we consider, prudentially, the sum of these two 
numbers. The estimate (|3.2|) is in substantial agreement with the MC estimate of 
Pc obtained using finite-size scaling techniques, /3c(A4 = 1.10) = 0.3750966(4). 
Similarly, for the (j)^ lattice model we obtain 

pciK = 1-90, Ag = 1) = 0.4269791(5), (3.4) 
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1.260 



1.255 



1.250 



n I I I I r 



\ I I I 

A spin-1 D=0.641 
* spin-1 /2 



1.245 



f 



1.240 



1.235 



1.230 



1.225 
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13 14 15 16 17 18 19 20 21 22 23 24 25 26 



Figure 2: Estimates of 7 as obtained by analyzing the HT series of x for the improved 
models and for the spin- 1/2 model, versus the order n of the series considered in the 
analysis. IA2's are considered. 



Table 3: Critical exponents obtained from the HT analysis. In parentheses we report 
the approximant error at A* or D*, in brackets the uncertainty due to the error on A* 
or D*, in braces the uncertainty due to the error on Pc- 











spin-1 


7 


X-scries 


1.23732(10)[12] 


1.23726(10) [22] 


1.23725(20) [10] 


a 


X-scrics (af) 


0.110(2) 


0.110(2) 


0.112(5) 


V 


^^-series 


0.6302(2) [1] 


0.6301 (3) [3] 


0.6300(2)[1] 


V 


^^-series (/3c-biased) 


0.63014(1){6}[9] 


0.63009(1){16}[16] 


0.63010(1){10}[9] 


r]v 


X,$^-series (CPRM) 


0.02294(3) [6] 


0.02291(2)[10] 


0.02294(8) [4] 
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le{K Ae = 1) = 1.23726(10) + 0.0055(A4 - 1.90), (3.5) 

and, using the MC result A* = 1.90(4), the estimate 7 = 1.23726(10) [22]; for the 
spin-1 model 

(3^{D = 0.641) = 0.3856717(10), (3.6) 

7e(L') = 1.23725(20) - 0.012(D - 0.641) (3.7) 

and therefore, using D* = 0.641(8), 7 = 1.23725(20)[10]. 

Our final estimate of 7 is obtained by combining the results of the three improved 
Hamiltonians: as estimate we take the weighted average of the three results, and as 
estimate of the uncertainty the smallest of the three errors. According to this rather 
subjective but reasonable procedure we obtain 

7 = 1.2373(2). (3.8) 

A direct estimate of the specific-heat exponent a is obtained from the singular be- 
havior of X at the antiferromagnetic critical point f3f = —^c, since R 



X = co + c,{P-pfy^' + ... (3.9) 

where 

^af = l-«. (3.10) 

FCIAk's provide rather precise estimates of ^af. The corresponding results for a are 
reported in the second line of Table ^ No error in brackets is reported since the 
dependence on A4, Z) is negligible. As final estimate we give 

a = 0.110(2). (3.11) 

The exponent u is obtained from the series of the second-moment correlation length 
^, since ~ {Pc — I3)~'^^ ■ Unbiased analyses of the 24th-order series of ^"^ / (3 provide 
the results reported in the third line of Table The corresponding estimates of (3^ 
are consistent with those derived from x, although less precise. For instance, for the 
04 model at A4 = 1.10 we found 13^ = 0.375098(2). 

In order to get a more precise estimate of u, we follow the procedure suggested in 
Ref. 0], i.e., we use the estimate of Pc obtained from x to bias the analysis of For 
this purpose we use lA's that have a singularity at a fixed value of f3c, or, in order 
to take into account the antiferromagnetic singularity, a pair of singularities at ±/?c; 
the two choices give equivalent results. This analysis provides the following effective 
exponents for the three classes of models. For A4 ~ A4 

i/e(A4) = 0.63014(1){6} + 0.0045(A4 - 1.10) (3.12) 

for the 0^ model, where the number in braces gives the variation of the estimate when 
(3c varies within one error bar; 

iy^(A4, Ae = 1) = 0.63009(1){16} + 0.004(A4 - 1.90) (3.13) 
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for the 0^ model; 

UeiD) = 0.63010(1){10} - O.Ollp - 0.641) (3.14) 

for the spin-1 modeL Then, using the MC estimates of A4, D*, one obtains the 
results reported in Table ^, where the error due to the uncertainty on A4 and D* 
is reported between brackets. They are perfect consistent with the results of the 
unbiased analysis, but more precise. Combining the results of Table |^ as we did for 
7, we obtain 

z/ = 0.63012(16). (3.15) 
Using the hyperscaling relation a = 2 — 3z^, we derive 

a = 0.1096(5), (3.16) 

which is fully consistent with, but more precise than, the direct estimate ( p.ll|) . 

Using the above- reported results for 7 and u and the scaling relation 7 = (2 — r/)//, 
we obtain rj = 0.0364(6), where the error is estimated by considering the errors on 7 
and u as independent, which is of course not true. We can obtain an estimate of rj with 
a smaller, yet reliable, error by applying the so-called critical-point renormalization 
method |^ to the series of x ci-nd This method provides an estimate for the 
combination rju. Proceeding as before, we obtain 

[riu]eiX4) = 0.02294(3) + 0.003(A4 - 1.10) (3.17) 

for the (j)'^ model, 

[w]e{K Ae = 1) = 0.02291(2) + 0.0025(A4 - 1.90) (3.18) 
for the 0^ model, and 

[w]e{D) = 0.02294(8) - 0.005(D - 0.641) (3.19) 

for the spin-1 model. We then obtain the results reported in Table |^, which lead to 
an estimate of 77 with a considerably smaller error: 

ri = 0.03639(15). (3.20) 

Then, by using the scaling relations we obtain 

6 = ^^ = 4.7893(8), (3.21) 
1 + 1] 

(3 = ^(1 + 77) = 0.32653(10), (3.22) 

where the error on (3 has been estimated by considering the errors of z/ and 1] as 
independent. 

Finally, we estimate the exponent A. For this purpose, we analyze the HT expan- 
sion of Vx that behaves like 

t^X = C-' {I + a^t^ + ...), (3.23) 
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for t = 1 — (3 / (3c ^ We consider the spin-1/2 model — here improved models are 
not useful since ^ — fix the exponent 7 to our best estimate, 7 = 1.2373, and use 
biased lA's that are singular at Pc = 0.22165459(10), which is the most precise MC 
estimate of the critical point We obtain 

A = 0.52(3), (3.24) 

where the error takes into account the uncertainty on j3c and 7. Correspondingly, we 
obtain u = A/u = 0.83(5). Consistent results are obtained from the analysis of the 
series of t'^'^^'^, fixing u and jSc- 



3.2 The ratio method 

In order to check the above-reported results, we consider the ratio method proposed by 
Zinn- Justin in Ref. (see also Ref. 0). Such a method has been recently employed 



in Refs. |llO|,[T7|| to analyze the 25th-order HT expansions of spin-S* models on the 
simple cubic and on the body-centered cubic lattice. 
According to this method, given a quantity 



S = 5^c„r ~ AsiPc - P)-'^ [1 + asiPc - PY + ...] 

n 

one considers the sequences 

/3W = f^!^zl^, exp 

where 



Sn + Sn-2 



n-2 



Sn-2) 



(3.25) 

(3.26) 
(3.27) 

_\n.{CnCn-i/ lll(Cn-lC„-5/c^_3)_ 

Asymptotically, the two sequences /3c"^ and approach /5c and with corrections 
of 0{l/n^^^) and 0{l/n'^) respectively. More precisely, if 

c„ ^ (Ao + A,n-^) (3.29) 

for n — >■ cxD, then 



[Sn — Sn-2) 
1 



+ 



pin) 



C 



in) 



C 



(3.30) 
(3.31) 



Note that, if only analytic corrections are present, i.e., e = 1, the convergence is faster 
with corrections of order and for Pc and C,: 

'AlC~2 7 _ .A 1 



Pc 

c 



1 - 



+ — (C- 1^ 

12^^ 





3,, 



Cn2 



(3.32) 

(3.33) 
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0.37512 



0.37511 - 



(n) 



0.37510 



0.37509 



0.37508 




Figure 3: Sequence (3^^^ for the 0^ model at A4 = 1.10 using the series for x- The 
dashed hues indicate the lA estimate of (3c- 



In Figs. ^ and ^ we show, for the 0^ model at A4 = 1.10 and for the spin-1 model 
at D = 0.641 respectively, the sequence /3c"^ obtained using S = x- The sequence 
clearly approaches the lA estimate. For the (p^ model the agreement is quite good and 
indeed /3c"^ differs from the lA estimate by 15 x 10"'' and 9 x 10"'' for n = 24, 25 
(note that the error on the lA estimate of /?c is 5 x 10^^). In principle, one could try 
to extrapolate the sequence /3c"^ to get a better estimate of (3c- For this purpose, we 
have tried to fit /J^"^ assuming a behavior of the form 



A 



a + bn 



(3.34) 



where a, b, and a are free parameters. If we interpolate /3c"'' for n = 21, 23, 25 with 
Eq. ( p.34| ), we obtain 



= 0.3750977(-5) + 3.0(+4) x 10" 



n N -6.6(-1.5) 



20 



(3.35) 



where the "errors" show the variation of the parameters between the interpolation 
with n = 21, 23, 25 and n = 19, 21, 23. Analogously, the even sequence n = 20, 22, 24 
gives 

/ n \ -6{-2) 

= 0.3750968(-25) + 4(+2) x 10"^ (^— j . (3.36) 

The extrapolated values are in perfect agreement with Eq. ( |3.2| ), but it is quite difficult 
to interpret the results for a. Indeed, in an improved model the leading corrections 
in the coefficients c„ are of order n"^^, n"^, with A2 ~ 1. The analytic term 
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0.38572 



0.38570 



0.38568 



(n) 



0.38566 



0.38564 - 



0.38562 



0.38560 




Figure 4: Sequence f3i^^ for the spin-l model at D 
The dashed hues indicate the lA estimate of (3c- 



0.641 using the series for x- 



gives a contribution of order n ^, while the nonanalytic one gives a correction of order 
j^-A2-i_ However, its amplitude is of order A2 — 1, and thus, since A2 ~ 1, it could 



be very small. The next correction terms are of order n 



n 



-l-A 



, and give rise 

to corrections of order n~^~^^, n~^~^. Inclusion of corrections with 2 < a < 5/2 
does not improve the fit |5^. Clearly, we are not yet sufficiently asymptotic to be 
able to extrapolate using the leading asymptotic behavior. At the values of n we are 
considering, several corrections are still important and apparently conspire to give a 
uniformly small correction. 

The same behavior is observed in the 0^ model, where both odd and even points 
extrapolate to 0.4269787, with effective exponent a ~ 12, 8. The agreement with the 
lA estimate ( p.4| ) is quite good. We finally analyze the spin-l results. Even points 
show again a very fast convergence with a ~ 9 and extrapolate to 0.3856662. Odd 



points instead are well fitted by assuming corrections of order n 01 n 



-5/2 



Fixing 



a = -2, we obtain 0.3856730, while for a = -5/2 we have 0.3856719. Again, the lA 
result ( |3.6|) is very well confirmed. 

For comparison, in Fig. |^ we plot the sequence /3c"^ for the spin-1/2 model versus 
1/n^/^ which should be approximately the leading correction. The higher-n results 
have apparently the predicted 0(n~^/^) behavior, and indeed an extrapolation with 
Eq. ( |3.34|) and a = 3/2 gives results that are close to the MC estimate of Pc- The 
odd (resp. even) points extrapolate to 0.22165686 (resp. 0.22165717): they are close 



to the MC estimate 0.22165459(10). However, it is hard to go beyond a relative 
precision of 10~^. 

In Fig. ^ we show the sequence 7*^"'^ as obtained from the series of x for the three 
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aooo 0.005 0.010 o.ois 0.020 0.025 



Figure 5: Sequence Pc for the spin-1/2 model using the series for x- The dashed 
hues indicate the MC estimate of (3c, while the dotted line corresponds to a n~^/^ 
extrapolation of the four points with n — 22, 23, 24, 25. 
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Figure 6: Sequences 7*^"^ for the (p* model at A4 = 1.10, the 0^ model at A4 = 1.90, 
the spin-1 model at D = 0.641, and the standard Ising model. The dashed hues 
indicate the lA estimate of 7. 
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Figure 7: Sequences [2z/](") for the (p'^ model at A4 = 1.10, the 0^ model at A4 = 1.90, 
the spin-1 model at D = 0.641, and the standard Ising model. The dashed lines 
indicate the lA estimate of 2i>. 



improved models and for the standard Ising model. The improved results clearly 
approach our best estimate 7 = 1.2373(2), the 0"^ and (p^ models from above and the 
spin-1 model from below. Note that the results are extremely flat and no extrapolation 
is needed. We also report the sequence 7*-"^ for the Ising model. If we extrapolate 
the results assuming a behavior of the form a + bn~^, with A = 0.52, we obtain 
7 = 1.23857, 1.23832, 1.23801 using pairs n = (21, 23), (22, 24), and (23, 25). Clearly, 
the estimates converge towards the lA estimate 7 = 1.2373(2). 

In Fig. 1^ we show the sequence [2;/]^"'' obtained from the series of Again, the 
improved models show a very good convergence to the lA estimate, in spite of the 
fact that the analysis is unbiased — the value of is not fixed. The Ising results are 
sensibly higher and steadily decreasing, reaching u ^ 0.638 for n = 25. Results that 
are closer to the lA estimate are obtained by an extrapolation. Assuming a behavior 
of the form a + bn~^, we obtain u = 0.6290 and 0.6284 from even and odd sequences 
respectively. Again, the agreement is satisfactory. 

In conclusion, this analysis based on the variant of the ratio method proposed by 
Zinn- Justin supports the lA estimate obtained in Sec. 

3.3 Matching the coefficients with their asymptotic form 

In the preceding section we have determined the critical exponents and i3c by gener- 
ating sequences that converge to the asymptotic value for 00. In this Section, 
following Ref. ||16|, we wish to perform a more straightforward analysis, both con- 
ceptually and practically. The idea is to generate sequences of estimates by fitting 
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the expansion coefficients with their asymptotic form. By adding a sufficiently large 
number of terms we can make the convergence as fast as possible, although of course 
the procedure becomes unstable if the number of terms included is too large com- 
pared to the number of available terms. In practice, one should include those terms 
that give rise to the maximal stability of the results. In some sense, the variant ratio 
method of the previous section corresponds to considering the leading singular be- 
havior and the ffist analytic correction — and also the leading nonanalytic term if we 
further extrapolate the sequence. 

On the cubic lattice, the large-order behavior is dictated by the singularities at 
±Pc- Indeed, given an observable S with expansion S = J^n^nP"", for n ^ oo the 
expansion coefficients behave like 



C-1 



Uo + 4 + ^ + ^ + 4 + ... 3.37 



Here, we have neglected all subleading exponents except the first one A, and in 
particular, the ffist subdominant A2. However, since A2 ~ 1 [0, for all practical 
purposes a term n~^^ cannot be distinguished from a purely analytic correction. 
Also, we do not write terms of order n~''^ since they cannot be distinguished from 
the analytic terms and corrections of order n~™~^. Note also the presence of the 
parity-dependent corrections with exponent ^af and the subleading corrections with 
exponent Aaf. For the susceptibility x 5 it is known |p3 that = 1— a. The argument 



can be generalized to all moments m2k and thus in all cases we predict 9i,f = 1 — a. 



We have tested this prediction for cf. Sec. |3.1| , m2, and 777,4. By analyzing the 
expansion of 7772 with biased lA/c's that have a pair of singularities in ±/3c, we obtain 
6'af = 0.884(12), while from the expansion of 7774 we obtain ^af = 0.90(9). These 
results are clearly compatible with the prediction Og^f = 1 — a = 0.8904(5). For the 
exponent Aaf nothing is known. We have analyzed the expansion using Aaf = 1/2 
and Aaf = 1. The results appear to be quite insensitive on either choice. For this 
reason, in the following we only report the results corresponding to purely analytic 
corrections, i.e., we set B^^^ = 0. 

Note that this method allows to determine the nonuniversal amplitudes Aq, A/^, 
. . ., and consequently the amplitudes appearing in the expansion of S* for /3 — > Pc- 
If 

S = Asi(3,-/3)-'^ [l + asil3,-(3f] (3.38) 

then 

As = r(C)Ao, (3.39) 

In the following, we shall perform two different analyses: (essentially) unbiased anal- 
yses in order to determine the exponents ( and (3c and biased analyses in which ( 
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and /3c are fixed. In all cases we fix the value of A, A = 0.52, and the exponent of 
the antiferromagnetic singularity. In the unbiased analyses, in order to have a linear 
problem, we consider lnc„ that behaves 

Inc. = -ln(/3,) n + iC- l)\nn + + ^ + ^ + + ^ + ■•• (3-41) 

As before, we have neglected terms that have exponents similar to those already 
present: for instance, terms 0{n~^^~'"^) or 0{n~^^'^~^^~™-) . In the expansion of the 
antiferromagnetic part we have assumed A^f = A, or A^f = 1. Note that if only 
analytic terms are present in Eq. ( |3.37|) , i.e., i?af = 0, then di is proportional to Ai 
and therefore it vanishes in improved models. 

We first analyze improved models and we verify that Ai ^ 0. For this purpose, 
we consider the susceptibility x and, for each improved Hamiltonian, we generate two 
sequences of amplitudes in the following way: 

(a) We choose two integers h, k and consider Eq. ( p.41|) keeping only Bq, . . ., bh-i 
in the ferromagnetic part and do, . . ., dk-i in the antiferromagnetic one. Then, 
we generate sequences /3c"\ 7*'"\ &o"^ • • ^i-D '^o"^ • • ^y solving the 

{h + k + 2) equations In Cn-m = Rn-m, m = 0, . . . , h + k + 1, where Rn is the 
right-hand side of Eq. ( p[TD . We use A = 0.52, 7 + ^af = 2.1277. 



(b) We choose two integers h, k and consider Eq. ( ^3.37] ) keeping only Aq, . . ., Ah-i 
in the ferromagnetic part and Bq, . . ., Bk_i in the antiferromagnetic one. We use 
A = 0.52, 7 = 1.2373, O^i = 0.8904, the lA estimate of fS^, and B^i = 0. Then, 
we generate sequences AI^\ . . ., bI"'\ . . ., -B["\, by solving the {h + k) 

equations c„_m = Rn-m, m = 0, . . . , h + k — 1, where Rn is the right-hand side 
ofEq. ( p37D . 

In both cases we vary h and k, trying to find the values that give the best stability of 
the exponents or of the leading amplitudes. In the unbiased analysis (a), the preferred 
choice is {h, k) = (4, 4), while for analysis (b) we use {h, k) = (3, 2). For these choices 
of the parameters, in Fig. we report the corresponding sequence of a^"'' = a^"\ 
obtained using Eq. (p.40|) . In the unbiased analysis (a), a^"^ clearly converges to zero 
for the improved Hamiltonians (f)^ and 0^, as expected. For the spin-1 model, it is not 
that clear, and presumably more orders are need to observe convincingly = 0. In 
the case of the biased analysis, a^"'' is very stable and small already for n > 15. For 
all Hamiltonians we observe \a^\ < 10~^. 

As a second check of consistency we have verified that our estimates of 
compatible with the quoted error bars on A4 and D*. For this purpose, using the 
analysis of type (b) reported above, we have computed for A4 ± AA4, where AA4 
is the quoted error bar — for the spin-1 model we are referring to D* ± AD. In all 
cases, we find Ictxl-^I ^ AA4)| > |a^(A4)| and that a^{Xl + AA4) and a^{Xl — AA4) 
have opposite sign. This confirms the correctness of our estimates of A4 and D*. Of 
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Figure 8: Amplitude of the leading scaling correction as obtained from several dif- 
ferent analyses of x for the standard spin-1/2 Ising model and for the three improved 
models. Details are explained in the text. 



course, since we use jSc and 7 obtained in the lA analysis, the above results represent 
only a check of consistency. Indeed: (i) we determine (3c and 7 by performing a lA 
analysis whose results should be reliable only if the models are improved (in some 
sense we assume weakly here ~ 0); (ii) using such values of jSc and 7, we estimate 
and find ~ 0. 

Once we have verified that Ai is very small and compatible with zero within the 
precision of the analysis, we have performed several analyses fixing Ai = and 61 = 0. 
At the same time, we have set di = 0, which corresponds to assuming Aaf = 1. We 
have determined the exponents by performing the analysis (a) reported above. In the 
case of the 0^ model for A4 = 1.10, this analysis gives 7 ~ 1.2374. Similarly, we obtain 
7 ^ 1.2375 for the 0® model at A4 = 1.90 and for the spin-1 model at D = 0.641. In 
Fig. ^we show the sequence j^^"^ for {h, k) = (5, 5) (since two coefficients vanish, we are 
considering four amplitudes in the ferromagnetic and antiferromagnetic expansion). 
We observe a very good agreement with the lA estimate 7 = 1.2373(2). It is difficult 
to estimate the uncertainty, since the results do not show a sufficiently robust stability 
with respect to the number {h, k) of coefficients used in the analysis. 

Finally, we report the estimates of the amplitudes obtained in the analysis of type 
(b) for the magnetic susceptibility: 

0^ : ^ 0.5246, A^^^ ^ 0.13, 5^^^ ^ -0.0351; 

06 : ^ 0.4601, A^^^ ^ 0.11, B^^^ ^ -0.0311; 
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Figure 9: Exponent 7 as obtained from the analysis (a) of x for the three improved 
models. Details are explained in the text. 



spin-1 : A'^^' ^ 0.5126, A'^' ^ 0.12, B^^> ^ -0.0359. 

Moreover, \A^^^\ < 10"^ for the 0^ and 0^ models, while ^3^^ ^ —0.02 for the spin-1 
model. Errors should be ±1 on the last reported digit and include the uncertainty 
on the n — >■ 00 extrapolation of the sequences and the variation of the estimates for 
h and k in the range h = 3 — 5 and k = 2 — 4. They do not take into account instead 
the variation of the estimates with 7 and jSc- Note that the estimate of A2 is purely 
phenomenological and in practice it should correspond to the sum of the amplitude 
of and of n~^^ (note that in improved models the amplitude of n~^^ vanishes). 

We have performed similar analyses for the spin-1/2 Ising model, in order to 
compute the nonuniversal amplitudes. We have performed: (a) an analysis of type 
(a) using {h, k) = (4,4); (b) an analysis of type (a) in which we have fixed f3c to its 
MC value using {h,k) = (4,3); (c) an analysis of type (b) using {h,k) = (3,2). The 
results for a^^ are reported in Fig. These analyses give perfectly consistent results 
and allow us to determine the amplitudes: 

a'i^'^ = 1.233 - 10(7- 1.2373) - 0.013(A- 0.52), (3.42) 
aS^) = -0.13 -0.7(A- 0.52) + 50(7 - 1.2373), (3.43) 
eJ^) = -0.073, (3.44) 

where we have explicitly written the dependence on the input parameters (when it 
turns out to be relevant). We have repeated the same analysis for the second moment 
1712. We obtain 

aJ,""^ = 1.301 - 10(7 + 2z/- 2.49754) - 0.07(A- 0.52), (3.45) 
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4"'^ = -0.73 - 4(A -0.52) + 55(7 + 2z/- 2.49754), (3.46) 
eJ'"^) = 0.06. (3.47) 



Using the above results and Eq. (p.40|) , one can determine the amphtudes and a^, 
associated with the O(t^) of the scahng corrections in the Wegner expansion of x and 
^ respectively, and evaluate their universal ratio. We obtain a^/a^ = 0.9(1), where 
the error takes also into account the uncertainty on the input parameters of the biased 
analysis. For comparison we mention the recent HT result a^/a^ = 0.76(6) Jl^, and 
the field theoretical estimate a^/cix = 0.68(2) ||3^ 



4 The critical equation of state 
4.1 Definitions 

The equation of state relates the magnetization M, the magnetic field H, and the 
reduced temperature t = {T — Tc)/Tc. In the neighborhood of the critical point t = 0, 
H = 0, it can be written in the scaling form 

H = B-'M'fix), (4.1) 
X = t{M/B)-^l^, (4.2) 

where B^. and B are the amplitudes of the magnetization on the critical isotherm and 
on the coexistence curve, 

M = B^H^'^ t = 0, (4.3) 

M = B{-tf H = 0, t<0. (4.4) 

Using these normalizations the coexistence curve corresponds to x = —1, and the 
universal function f{x) satisfies /(—I) = 0, /(O) = 1. Griffiths' analyticity implies 
that f{x) is regular everywhere for x > — 1. It has a regular expansion in powers of 

X, 

oo 

/(x) = l + 5^^V, (4.5) 

n=l 

and a large- a; expansion of the form 

oo 

/(x)=x^5^/rx-2"^. (4.6) 



n=0 



At the coexistence curve, i.e., for x —1, f{x) has at most an essential singularity 
||58| . It can be asymptotically expanded as 

oo 
n=l 
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It is useful to rewrite the equation of state in terms of a variable proportional to 
Mt~^, although in this case we must distinguish between t > and t < 0. For t > 
we define 



H 



z = 



ct) '"'^'^ 



ci 



1/2 



Mt- 



(4.8) 



while for t < we set 



(4.9) 



The constants and are the amplitudes appearing in the critical behavior of 
the zero- momentum connected n-point correlation functions Xn'- 

The susceptibility x corresponds to X2 and we simply write = C^. 



With the chosen normalizations |4T|,^,^ 
F(z) = 



1 3 1 

H — z + > —roj z 

6 ^(2j-l)! 



2i-i 



i=3 

oo 



i=3 '• 



The functions F{z^ and $(«) are related to f{x). Indeed, 



and 



z''F{z) = F^fix 



(4.11) 
(4.12) 

(4.13) 
(4.14) 

(4.15) 

To compare with experimental data, it is useful to determine the magnetization as a 
function of tH~^/^^ . Therefore, we define 



Z = ZqX ^, 



U = —X 



The constant is defined by the large-z behavior of F{z) 



F{z)=z'J2^> 



while 



k=0 



1/2 



B. 



E{y) = B;'MH-^/^ = f{x)~^'\ 

y ^ (5/5,)V/3^^-i/W = a;/(a;)-i/(^^\ 



(4.17) 
(4.18) 
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Finally, we shall also determine the scaling behavior of the susceptibility, by defining 



4.2 Small-magnetization behavior 

In this section we determine the first few coefficients r2j appearing in the expansion 
of the scaling function F{z), cf. Eq. ( [4 .81) . We shall also compute the four-point 
renormalized coupling constant g^, which, although not related to the equation of 
state, is relevant for the field-theoretical approach and will be used to determine 
amplitude ratios involving the second-moment correlation length. 

In order to estimate the critical limit of and of we first determine their HT 
expansions using the corresponding results for X2j and m2 

9. = (4.20) 
re = 10-^, (4.21) 

xi 

2 

rs = 280 - 56^ + (4.22) 

xi xi 



The corresponding series |59] have been analyzed following closely the procedure 



presented in App. B.3 of Ref. [^. We use biased lAl's with a singularity at jSc or a 
pair of singularities at ±Pc, where Pc is obtained from the analysis of the susceptibility. 
Around Pc, lAl's behave like 



lAl^f{(3){l-P/P,)^ + g{P), (4.23) 

where /(/3) and g{P) are regular at Pc, provided ( is not a negative integer. In 
particular 

(see Eq. ( pTTP for the definition of the above quantities). In the case we are considering, 
C is positive and, therefore, g{Pc) provides the desired estimate. 

In Table ^ (first line) we report the estimates of g^ obtained for the three improved 
Hamiltonians. The error in parentheses is related to the spread of the approximants 
and the second one in brackets to the uncertainty on Xl,D*. The error induced by 
the uncertainty on Pc is negligible. The results are perfectly consistent. Our final 
estimate is 

^4 = 23.56(2). (4.25) 

The result for the exponent ( in Eq. (|4.23| ) is C = 1.3(3), which is consistent with 
our expectation for improved models, i.e., C = ^2 ~ 2 A and A ^ 0.5. For com- 
parison, the same analysis applied to the standard Ising model gives g^ = 23.5(5) 
and ( = 0.6(3), in agreement with the fact that in this case ( = A. Notice that 
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Table 


4: Results 


for g4, rg and 


^^8- 








spin-1 


final estimates 


94 


23.559(8)[11] 


23.554(8)[20" 


23.560(20) [5] 


23.56(2) 


re 


2.057(4) [1] 


2.056(4)[2] 


2.052(8)[2] 


2.056(5) 


rs 


2.29(9)[3] 


2.31(5)[5] 


2.37(7)[3] 


2.3(1) 



the small difference with the estimate of reported in Ref. [|T9| is essentially due to 
the different analysis employed here, which is better justified due to the nonanalytic 
behavior at j3c predicted by renormalization group fHlI]. With respect to standard 
Fade approximants, biased lAl's require more terms of the series to give reasonable 
results, but they are less subject to systematic errors since they allow for confluent 
nonanalytic corrections at Pc- Biased lAl's give (74 = 23.54(4) when applied to the 
17th-order series of Ref. [|l^ . 



Results for re, rg are obtained using the same method and are reported in Table ^ 
We finally recall that a rough estimate of rio was obtained in Ref. |jl9| from the 
analysis of its 15th-order series, obtaining rio = —13(4). A review of the available 
results for these quantities can be found in Ref. [|ll. 



4.3 Parametric representations of the equation of state 

In this section we shall determine the equation of state using parametric representa- 
tions, improving the results of Refs. |T9|,^. This method has also been applied in 
two dimensions and to the three-dimensional XY p5|,p9| and Heisenberg 
universality classes. 

In order to obtain approximate expressions for the equation of state, we parametrize 
the thermodynamic variables in terms of two parameters R and 6, implementing all 
expected scaling and analytic properties. Explicitly, we write ||63|-|65| 

M = moR'^e, 

t = R{l-e^), 
H = hoR^^h{e), (4.26) 

where ho and mo are normalization constants. The function h{6) is odd and normal- 
ized so that h{6) = 6 + 0{6^). The smallest positive zero of h{6), which should satisfy 
6^0 > 1, corresponds to the coexistence curve, i.e., to T < Tc and H 0. We mention 
that alternative versions of the parametric representations have been considered in 
Ref. H. 

in terms of 6. 



It is easy to express the scaling functions introduced in Sec. ff] 
The scaling function f{x) is obtained from 



X 



1 - ^2 fe, 
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/(^) - r'M, (4.27) 

while F{z) is obtained by 

z = p9 {1-9')-', 

F{z{e)) = p{i-eyh{e), (4.28) 

where p can be related to mg, Hq, and 6*4^ using Eqs. ( [4 .81 ) and ( [4.2(j| ). 

It is important to note that Eq. ( [4.26| ) and the normalization condition h{9) ~ 9 
for ^ do not completely fix the function h{9). Indeed, one can rewrite the relation 
between x and 9 in the form 

= h{l) /o°°(l - 9'y9^-\ (4.29) 

Thus, given /(x), the value of h{l) can be arbitrarily chosen to completely fix h{9). 
In the expression ( [4.28|) we can fix this arbitrariness by choosing arbitrarily the pa- 
rameter p [|T],[T|,^,|T[]. 

As suggested by arguments based on the e-expansion |]T9|,^, we approximate h{9) 
with polynomials, i.e., we set 

k 

h{9)=9 + Y,h2n+i9''^^'. (4.30) 

n=l 

This choice is further supported by the effectiveness of its simplest version with k = 
1, which is the so-called linear model. If we require the approximate parametric 
representation to give the correct {k — 1) universal ratios rg, rg, . . ., r2k+2, we obtain 

n 

h2n+l = Cnn.6™(/13 + 7)™7t/^^^, (4.31) 

^ (2m +1)! 

m=0 ^ ' 

where 

^ n—m 

Cnm = j TT TT(2/3m-7 + A;-l), (4.32) 

(n-m)! 

and we have set r2 = r4 = 1. Moreover, by requiring that F{z) = z + + we 
obtain the relation 

p2 = 6(/i3 + 7)- (4.33) 

In the exact parametric representation, the coefficient can be chosen arbitrarily. Of 
course, this is no longer true when we use our truncated function h{9), and the related 
approximate function /ipprox(a;, h^) depends on h^. We must thus fix a particular value 
for this parameter. Here we use a variational approach, requiring the approximate 
function /ipprox(x, /is) to have the smallest possible dependence on h^. Thus, we set 
hs = hs^k, where h^^k is a solution of the global stationarity condition 



dh: 



= (4.34) 

'13=^3, fe 
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Table 5: Polynomial approximations of h{6) using the global stationarity condition 
for various values of the parameter k. The reported expressions are obtained by using 
the central values of the input parameters. The last column shows the corrections to 
the simple linear model /iiin(6', ^o) = ^(1 — ^^/^o)- 



k 




h{d)/0 






hi9)/hun{o,eo) 


1 


1 


- 0.7347326'^ 


1.36104 


1 




2 


1 


- 0.731630612 + 0.00909004 


1.39085 


1 


- 0.01264296*2 


3 


1 


- 0.7367436(2 + 0.0089046l4 - 0.0004726*^ 


1.37861 


1 


- 0.01137756*2 + 0.00065116*4 



for all X. Equivalently one may require that, for any universal ratio R that can be 
obtained from the equation of state, its approximate expression i?ipprox obtained by 
using the parametric representation satisfies 



dR 



(fc) 

approx 



dh?, 



0. 



(4.35) 



The existence of such a value of is a nontrivial mathematical fact. The stationary 



value of /13 is the solution of the algebraic equation |T9| 



2(2/3-l)(/i3 + 7)^-27 + 2A; 



h 



2k+l 



0. 



(4.36) 



For /c = 1, the so-called linear model, Eq. ( [4. 361) gives 

7(1 - 2(3) 



7-2/3 



(4.37) 



which is the optimal value of considered in Ref. [^. Thus, the optimal (sometimes 
called restricted) linear model represents the first approximation of our scheme. 



4.4 Results 

Following Ref. |jT9|, we apply the variational method by using the HT results for 
7 = 1.2373(2), u = 0.63012(16), rg = 2.056(5), rg = 2.3(1), and r^o = -13(4) as input 
parameters of the approximation scheme. This provides different approximations with 
k = 1,2,3,4. In Table ^ we report the polynomials h{6) for k = 1,2,3, that are 
obtained in the variational approach for the central values of the input parameters. 
The fast decrease of the coefficients of the higher-order terms in h{6) gives further 
support to the effectiveness of the approximation scheme. We do not report h{6) 
for k = 4, since it requires rio and its available estimate is rather imprecise. Using 
the results reported in Table ^ and Eqs. ( [4.27] ), ( |4.28| ), and ( f4.14| ), one may easily 



27 




compute the corresponding approximations for the scahng functions f{x), F{z), and 
The results show a good convergence with increasing k. Actually, the results 
for k = 2,3,4 are already consistent within the errors induced by the uncertainty on 
the input parameters, indicating that the systematic error due to the truncation is at 
most of the same order of the error induced by the input data. In Figs. 0, [IT], and 



12] we show respectively the scaling functions as obtained from h{9) for k = 1,2, 3. 
In Table ^ we report results concerning the behavior of the scaling function f{x) 



F{z) and ioi H = and on the critical isotherm, cf. Eqs. (4^), (O), ( |4.7| 



( pnp , ( CT ), (|4l5| ). Note that the results for k = 1,2,3 oscillate and that the 



uncertainty due to the input parameters on the k = 3 results is approximately the 
same as the difference between the estimates with k = 2 and k = 3. Therefore, 
it is reasonable to consider the k = 3 truncation as the best approximation of the 
method using the available input parameters and to use the corresponding errors as 
final uncertainties. 



In Fig. we give the behavior of the magnetization as a function of t and 
H, reporting the scaling function E{y). The behavior of the susceptibility can be 
obtained from the scaling function D{y). The function D{y) has a maximum for 
1/max = 1.980(4), corresponding to the so-called crossover or pseudocritical line (see 
Sec. 1^). In order to simplify possible comparisons, it may be convenient to consider 
the rescaled function 

D{y) 
C{yR) = jTf V' 
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Figure 11: The scaling function F{z). We also show the plot of the small- 2; expansion 
(dotted line), i.e., F{z) ^ z + ^z^ + -^rez^ for z ^ 0. 




1.0 1.5 2.0 2.5 

u 



Figure 12: The scaling function We plot also the asymptotic behavior of 

at the coexistence curve (dotted line), i.e., ^{u) ^ {u — l) + ^V3{u — 1)^ + ^V4{u — 1)^ 

for M — > 1. 
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Table 6: Expansion coefficients for the scaling equation of state obtained by the 
variational approach. See text for definitions. Numbers marked with an asterisk are 
inputs, not predictions. 





k = 1 


fc = 2 


A: = 3 


/fc = 4 




1.3610(8) 


1.390(2) 


1.38(2) 


1.34(6) 


o 

r 


1.7365(8) 


1 741fl) 


1.733(10) 


1.69(6) 


' o 


1.938(3) 


*2. 056(5) 


*2. 056(5) 


*2. 056(5) 


rs 


2.50(2) 


2.39(3) 


*2.3(1) 


*2.3(1) 


rio 


-12.59(2) 


-12.08(5) 


-10.6(1.8) 


*-13(4) 


Zpoo 


0.03277(8) 


0.03388(11) 


0.03382(15) 


0.0338(2) 


zo 


2.8254(7) 


2.792(2) 


2.794(3) 


2.798(8) 


f? 


1.05041(7) 


1.0532(2) 


1.0527(7) 


1.051(2) 




0.04298(6) 


0.04494(13) 


0.0446(4) 


0.0439(13) 




-0.02474(4) 


-0.02595(8) 


-0.0254(7) 


-0.023(4) 


£00 

Jo 


0.5960(4) 


0.6031(7) 


0.6024(15) 


0.601(4) 


rcocx 

Jl 


0.93912(9) 


0.9347(3) 


0.9357(11) 


0.938(4) 




6.013(4) 


6.062(4) 


6.050(13) 


6.02(5) 




16.32(3) 


16.10(4) 


16.17(10) 


16.4(3) 



Vr 



(4.38) 



which is such that the maximum corresponds to = 1 and satisfies C(l) = 1. In 
Fig. |14| we plot the scaling function C{yji), as obtained from the k = 1,2,3 approxi- 



mate parametric representations. 

In experimental work on magnetic systems, it is customary to report [^,^ h/m = 
H\t\~'^/M versus = M'^\t\~'^^. Such a function can be easily obtained from our 
approximations for f{x), since = B'^\x\~'^^ and 



— = k\x\ ^f{x), 
m 



where the constant k can be written as 



k = B-^B^'/^ = 



(4.39) 



(4.40) 



where = C^B^^^ /B^ is a universal constant, see Sec. ^ A plot of w? / B"^ versus 
C^h/m for the two cases t > and t < is reported in Fig. |l^. 

It is interesting to observe that in a neighborhood of the critical isotherm the 
equation of state can be written in the Arrott-Noakes form 



H 

M 



1/7 



at + bM^'^, 



(4.41) 
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Figure 15: Plot of w? jB'^ versus C^h/m. 



where a and h are numerical constants. Indeed, using the results of Table ^ for = 3, 
we obtain 



H 



M 



1/7 



_ fc-V7 = _ + 0.851 t - 0.050 — 



M 



B 



-1//3 



0.008 ( ^"j 



-2/3 



(4.42) 

Thus, corrections to Eq. ( [4.41| ) are small and thus this expression has a quite wide 
range of validity. 



5 Universal amplitude ratios 

From the critical equation of state one may derive estimates of several universal am- 
plitude ratios. They are expressed in terms of the amplitudes of the magnetization, cf. 
Eqs. ( |4.3D and (|4.4| ), of the magnetic susceptibihty and n-point correlation functions, 
cf. Eq. (|4.10| ), of the specific heat 

CH = A^\t\-'^, (5.1) 
of the second-moment correlation length 

e = /^lr^ (5.2) 

and of the true (on-shell) correlation length, describing the large distance behavior 
of the two-point function, 

^gap = /gapW • (5-3) 
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^ u 


= A+ /A- 


Uo = 

^2 — 




Ua 


= Ct /C7 


Rt = 


^4 ^ 1 1 


R+ 

c 


= aA+C+/B^ 


i?- = 

c 


aA'C-/B^ 


R7 

± 




R. = 


C+B^-^HBrY 






- 


^4 1 \^ 1 ' 


Hi 


^4 / LV J \J 1 \ 




^ / L Vj / J 


Ue 


— J 13 


1 ■ 


= f+ /f" 
— J gap/ J gap 




= aA+(f+)^ 


0" = 


aA-(f-Y 








~ t/g£ip / if' 




^ gap / f 


Ql^ 


fc irc 
J gap / J 


Qc 


^ B\f+f/C+ 


Q2 = 


(fyf+f-vc+/c. 


p 




Pc = 




Rp 


= C+/Cp 







Table 7: Amplitude-ratio definitions. 



One can also define amplitudes along the critical isotherm, e.g. 

^ = C'\H\-^/f^^, (5.4) 

e = mr/^', (5.5) 

^gap = fgap\H\ . (5-6) 

We also consider the crossover (or pseudocritical) line tmax(-f^), that is defined as 
the reduced temperature for which the magnetic susceptibility has a maximum at H 
fixed. Renormalization group predicts 

t^,^{H)=TpH^/^^^^\ (5.7) 
x(Wx,i^) = Cpt-L. (5.8) 

We consider several universal amplitude ratios. They are defined in Table ^. 

In Table |]we report the universal amplitude ratios, as derived by the approximate 
polynomial representations of the equation of state for k = 1,2,3,4. The reported 
errors are only due to the uncertainty of the input parameters and do not include 
the systematic error of the procedure, which may be determined by comparing the 
results of the various approximations. In Table |^ we also show results for Zraa.x, a^max 
and Umax which are the values of the scaling variable z, x and y {y was defined in 
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Table 8: Universal amplitude ratios obtained by taking different approximations of 
the parametric function h{6). 





k = 1 


k = 2 


k = 3 


fc = 4 


Uo 


0.5231(11) 


0.533(2) 


0.5319(25) 


0.529(6) 


U2 


4.826(7) 


4.745(10) 


4.758(19) 


4.78(5) 


Ui 


-9.73(3) 


-8.85(6) 


-9.0(2) 


-9.3(5) 


Rt 


0.05545(7) 


0.0570(1) 


0.0567(3) 


0.0562(11) 


Rc 


0.021967(11) 


0.02253(3) 


0.02242(12) 


0.0222(4) 


Rt 


7.983(4) 


7.794(8) 


7.81(2) 


7.83(4) 


R4 


92.15(13) 


94.13(13) 


93.6(6) 


92(2) 


Rx 


1.6779(11) 


1.658(2) 


1.660(4) 


1.665(10) 


U2Rt 


38.52(5) 


36.98(10) 


37.1(2) 


37.4(6) 




0.4427(7) 


0.4444(7) 


0.443(2) 


0.440(6) 


Pra 


1.25203(6) 


1.2493(2) 


1.2498(6) 


1.251(2) 


Pc 


0.3831(3) 


0.3938(5) 


0.3933(7) 


0.3930(11) 


Rp 


1.9789(3) 


1.9658(6) 


1.9665(10) 


1.9671(16) 




1.2443(4) 


1.2317(5) 


1.2322(8) 


1.2326(12) 




12.32(3) 


12.26(3) 


12.27(4) 


12.31(8) 


^max 


1.990(1) 


1.977(2) 


1.980(4) 


1.984(9) 


£'(Wniax) 


0.36179(4) 


0.36277(7) 


0.36268(14) 


0.3626(3) 



Eq. ( |4.18| )) associated with the crossover line. As already mentioned in Sec. |4.4] , we 
consider the k = 3 results as our best estimates, and report them in Table |I|. 

Estimates of universal ratios of amplitudes involving correlation-length ampli- 
tudes, such as Q~^, R'^, and Qc, can be obtained using the HT estimate of g^. For 
instance = R^R^ /g^. Other universal ratios can be obtained by supplementing 
the above results with the estimates of and (see Table ^, obtained by an anal- 
ysis of the corresponding low-temperature expansions |14,0|, and the HT estimate 



of (see Sec. H). Moreover, using approximate parametric representations of the 
correlation length, see Refs. [|l],|l9| for details, one may also estimate the universal 
ratios and Q2 defined in Table 0. 

6 Low-momentum behavior of the structure factor 

In this section we update the determination of the first few coefficients that parametrize 
the low-momentum expansion of the scaling two-point function in the HT phase 

00 

g{y)^^ = l+y + Y,(^,y\ (6.1) 

where y = k'^^^. 

The coefficients Cj can be related to the critical limit of appropriate dimensionless 
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Table 9: Estimates of the coefficients q, z = 2, 3, 4, of the low-momentum expansion 
of the structure factor. 











spin- 1 


final estimates 


C2 


-0.390(7) X 10-^ 


-0.390(6) 


X 10-'' 


-0.389(12) X 10-^ 


-0.390(6) X 10-^ 


C3 


0.882(8) X 10-5 


0.882(6) X 


10-5 


0.88(4) X 10-5 


0.882(6) X 10-5 


C4 


-0.4(1) X 10-6 


-0.4(1) X 


10-6 


-0.4(1) X 10-6 


-0.4(1) X 10-6 



ratios of spherical moments m2j. See Ref. for details. We have estimated the 
first few coefficients Cj from the corresponding series derived from the 25th-order 
expansions of m2j, using the analysis described in Sec. |4.2| . The results for the three 
improved models and our final estimates are reported in Table |^. Other interesting 
quantities are 



Sm = ^ga.p, 
Sz = XM^Z, 



:ap) 



(6.2) 
(6.3) 



where Mgap (the mass gap of the theory) and Zgap determine the long-distance be- 
havior of the two-point function: 



G{x) 



\x\ 



4:Tt\x\ 



(6.4) 



As discussed in Refs. |T9|,0, one may estimate Sm and Sz from C2, C3, and C4. Indeed, 
we have 



Sm = 1 + C2 - C3 + C4 + 2C2 + ... 
Sz = l-2c2 + 3c3-4c4-2c2 + 



(6.5) 
(6.6) 



where the ellipses indicate contributions that are negligible with respect to C4. There- 
fore, one finds Sm = 0.999601(6) and Sz = 1.000810(13). From the result for Sm, 
one obtains Q+ = f+p/f+ = 1.000200(3). 

A more detailed analysis of the behavior of the structure factor for all momenta 
can be found in Ref. . 
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